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Superconducting trapped-field magnets: Temperature and field distributions during 

pulsed field activation 
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We calculate the temperature and magnetic field distributions in a bulk superconductor during the 
process of pulsed-field magnetic activation. The calculations are based on the heat diffusion equation 
with account of the heat produced by flux motion, and the critical state model with temperature 
dependent critical current density. For a given activation time the total amount of trapped flux <1> 
is maximum for an optimal value Bopt of the maximal applied field. We analyze how Bopt and <& 
depend on the material parameters and the field ramp rate. 
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I. INTRODUCTION 

The recent progress in fabrication of large sized high- 
temperature superconductors with high critical current 
density makes them extremely promising for use as per- 
manent magnets. Trapped fields exceeding 12 T have 
already been reported for a YBa2Cu307_5 magnet at 
22 K.EI To magnetically activate the superconductor one 
most often uses a pulsed field magnetization (PFM). 
Whereas the PFM activation method is the most con- 
venient from practical point of view, even higher J 
have been trapped by quasi-static field ramping.E 
This shortcoming of PFM cannot be overcome by simply 
increasing the maximal applied field Bm- In fact, it is 
found experimentally that the trapped flux reaches max- 
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This behavior is believed to re- 
sult from heating produced by flux motion, which leads 
to strong temperature rise in superconductor during the 
activation process. There have already been suggested 
practical ways to improve tho-siUiation, in particular, by 
applying multiple field pulses .qeHj Even better results can 
be expected if the development is accompanied by a mod- 
eling of how the heat actually dissipates and redistributes 
in a superconductor during the PFM. Surprisingly, al- 
most no efforts have been done in this direction, and this 
lack of insight motivates the present work. 

Our analysis is based on the heat diffusion equation, 
taking into account the time and position dependent heat 
dissipation due to flux motion. As a result, the tem- 
perature and magnetic field distributions in the super- 
conductor during all the stages of the PFM process are 
calculated. We determine the optimal applied field Bopt 
corresponding to the maximum trapped flux <i>(-Bopt) and 
analyze how Bopt and <l>(_Bopt) depend on the parameters 
of superconductor. 



V X B = ^oj , 

V X E = -dB/dt. 



(1) 
(2) 



These equations should be supplemented by a relation- 
ship between j and the flelds E and B, which depends on 
the superconductor material, as well as on temperature T. 
When a superconductor is subjected to a non-stationary 
external magnetic fleld, Ba{t), a heat per unit volume is 
produced with the rate 



= Ej. 



(3) 



The heat release creates a temperature rise which is given 
by the thermal diffusion equation 



CdT/dt- kV^T = W . 



(4) 



Here k is the thermal conductivity, and C the heat ca- 
pacity per volume. 

For simplicity we will make calculations for a supercon- 
ductor shaped as a slab, while we expect that all qual- 
itative results are valid also for other geometries. The 
superconductor occupies space \x\ < w, and satisfles the 
boundary condition T(|a;| = w) = Tq, which assumes 
ideal thermal contact with the surroundings at the tem- 
perature Tq. The solution for a uniform initial tempera- 
ture can be expressed as 



dx' I dt' G{x,t-x' ,t') 

-w Jo 



n W{x',t') 



C 



(5) 



where G{x, t; x' , t') is the Green's function due to a unit 
instantaneous plane source at x' at t' . The Green's 
function satisfying the boundary conditions G{x — 
±w, t; t', x') = is given byll3 



II. EQUATIONS AND BASIC APPROACH 

The distributions of magnetic field, B, electric field E 
and current density j inside a sample are determined by 
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To calculate heat release rate W{x,t) we use the crit- 
ical state model, according to which the magnitude of 
the current density in flux-penetrated regions of a super- 
conductor equals to some critical value, jc- The critical 
current density depends in general on both the local field 
B and local temperature T, thus the magnetic field pro- 
file is determined by Eq. (|l|) with |j| = jdB, T). In order 
to obtain analytical results, we choose the Bean model, 
i. e. assume jc to be B independent. To account for the T 
dependence of jc the following iterative procedure is used. 
First, jc is taken T-independent and the time evolution 
of the profiles B{x, t) and W{x, t) is calculated. The W is 
substituted in Eq. (^) to determine the temperature pro- 
file T{x,t). It is then used to recalculate the jc[T{x,t)] 
which subsequently gives corrected magnetic field profiles 
according to Eq. ffl). Fortunately, it turns out that even 
the 1st iteration gives correct results within very good 
accuracy for realistic parameters. This is demonstrated 
below by a self-consistent numerical solution of the above 
set of equations. 

Let us consider a PFM where the external field is ap- 
plied as a triangular pulse. 

Bait) = R {tM ~ \t - tM\) , (7) 

where i? is a constant ramp rate, i. e. the field increases 
to the maximum value Bm during <t < Im ~ B^/R^ 
and then decreases to zero during Im < t < 2tM- When 
the temperature dependence of the critical current is ne- 
glected the field profile is given by the conventional Bean 
model for a zero-field-cooled superconductor, i.e., when 
the external field increases, the magnetic flux occupies 
the region Xf){t) < \x\ < w, where 

xo{t) - vt, V = R/iiojco , jco = jc(Tb) . (8) 

The heat release in the penetrated region is obtained from 
Eqs. d) and I): 

W{x,t) = [\x\-xo{t)]R3co (9) 

while W{x,t) = at |a;| < xoit). After the sample be- 
comes fully penetrated, i. e. at t > w/v, one has 

W{x,t)^\x\Rjco. (10) 

Similarly, for decreasing applied field we obtain 

W(x,t)^[\x\-xi{t)]RjcQ, xi{t)<\x\<w, (11) 

where 

xi{t)=w-v{t-tM)/2 (12) 

is the position of maximum flux density. In the region 
a; I < xi{t) there is no flux motion, and therefore no 
heat release. For t > tM + ^w/v the field is decreas- 
ing throughout the sample, and the dissipation is again 
given by Eq. (^0|). From this set of W{x,t) one finds the 
temperature distributions T{x, t) at all stages of the pro- 
cess. The expressions are listed in the Appendix. Finally, 



we determine the refined B{x,t)) by assuming a linear T 
dependence of the critical current, 

Tc — T{x, t) 
jc(x,t) = jcQ Tf; — (13) 

where Tc is the critical temperature.0 

The present thermo-magnetic problem is characterized 
by only two dimensionless parameters. The first is the 
ratio of Im to the thermal diffusion time CP /k, 

T = tMK/Cf = Bmk/RCw^ . (14) 

If T ^ 1 the heat diffusion can be neglected, whereas 
for T 3> 1 the heat escapes the sample so fast that the 
temperature increase is negligible. The second parameter 
is 

Bp ( fiQw'^jco djcV^^ , . 

^=Bir[-^of) • 

Here Bp = fJ-ojco'w is the full penetrjation field, and Bfj is 
the threshold field for a flux jump.ll3ll3 Consideration of 
flux jumps - macroscopic flux avalanches accompanied by 
pronounced heating - is beyond the scope of the present 
study. In practice, one wants to avoid flux jumps, which 
ruin the magnetization process, and can even damage the 
material. Therefore we limit ourselves to the parameter 
range where a < 1, and where the above iteration proce- 
dure is applicable. 

III. RESULTS AND DISCUSSION 

We present results of the calculations assuming a sam- 
ple of size 2u' = 1 cm and characterized by jco = 4.4 • 10^ 
A/m2 (giving Bp ^ 2.8 T), C = 0.88 x 10^ J/m^K, 
and K = 6 J/msK, all of which are typical values for 
YBa2Cu307_5 crystal with = 92 K at To = 77 K. 
From Eq. (^3|) one finds that the above set of parameters 
corresponds to a = 0.48 < 1. 

Shown in Figs. |l| and ^ are the temperature and field 
distributions after completed PFM, i.e., in the remanent 
state at t — 2tM- The figures display the result for the 
cases Bm = 2Bp and _Bj\/ = l.QBp, and the profiles are 
plotted for four different ramp rates of the applied field, 
corresponding to r ranging from 1 to 0.001. At high rates 
the temperature profiles have pronounced maxima near 
the sample surface. The reason is obvious - near the 
surface the fiux motion is most intense and the heat re- 
lease is maximal, see Eqs. (|)-(0. For a high ramp rate 
(r <C 1), the heat remains mostly in the regions where 
it was released. If the rate is low (t > 1) the heat has 
time to diffuse both to the center and the surface. As a 
result, T{x,t) becomes more uniform, but on the average 
the temperature rise is smaller than for high rates. 

As seen from the lower panels the temperature rise 
strongly affects the profiles of the remanent field. In- 
creased temperatures give a lower critical current, and 
hence, less steep slopes \dB/dx\ — jiojciT). For the high- 
est ramp rate the profiles are also most non-linear, and 



3 




-1.0 -0.5 0.0 0.5 1.0 

X / w 



FIG. 1: Temperature and flux density profiles in the remanent 
state after application of magnetic fleld Bai ~ 2Bp. Four 
curves are calculated for different durations of the field pulse: 
2tM ~ 7.3 s (solid), 0.73 s (dashed), 73 ms (dotted) and 7.3 ms 
(dash-dotted), which correspond to r = 1, 0.1, 0.01 and 0.001, 
respectively. 

the total amount of trapped flux is lowest. Hence, we 
reproduce the expected result that very low ramp rates 
give maximum trapped flux. 

In practice, the duration of the fleld pulse is limited. 
The key point is then to choose the optimum applied field 
Bm so that in the remanent state the trapped flux be- 
comes maximum. If Bm is too small the flux penetration 
is also small, whereas too large Bm gives an excessive 
heating and little flux becomes trapped by the sample. 
This is illustrated in Fig. H, where the total amount of 
trapped flux $ = / B{x)dx per unit length of the slab is 
plotted as a function of Bm- The solid curves show the 
conventional Bean- model result, 

$ = wBlj/{2Bp), Bm < Bp, 

= w [2Bm -Bp- Blj/{2Bp)] , Bp < Bm < 2Bp, 
= wBp, Bm > 2Bp, (16) 

obtained for a slab superconductor assuming linear field 
profiles. This result is applicable when the heating is 
negligible, i.e., for very low ramp rates or large C or k, 
i.e., when r — > oo or a — > 0. The presence of heating leads 
to reduction of the trapped flux, and a peak in the ^(Bm) 
curve appears at some Bm between Bp and 2Bp. The two 
panels of Fig. ^ allow us to trace the effect of changing 
the parameters r and a. Plotted as symbols in Fig. y 
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FIG. 2: The same curves as in Fig. y but for the maximum 
applied field Bm ~ l.GBp. 



are the exact results obtained by numerical simulations, 
whereas the lines represent the analytical solution. For 
the most important range of Bm where $ has the peak, 
the agreement is excellent, and only a slight deviation 
appears for large Bm- We conclude therefore that the 
model given by Eqs.(|), (|), allows to determine 

the Bm producing the maximum trapped flux. 

Let us consider 3 characteristic points, A, B and C on 
the $(Bm) where the peak is most pronounced, and an- 
alyze the corresponding B and T distributions. These 
remanent profiles are shown in Fig. ^. For small Bm 
(A) too little fiux penetrated the superconductor, and 
the B{x) profile has a large dip in the center region. For 
large Bm (C), the B{x) profile has the "right" triangular 
shape, however, its slopes are not maximally steep due 
to heating. In the state of maximum trapped flux (B) 
the small dip in the center is compensated by having rel- 
atively large slopes in the overall peak. Interestingly, we 
find that the optimum case always has a small minimum 
in the flux density at the center of the superconductor. 

If the maximum applied field is very large, the heat is 
most of the time released uniformly throughout the super- 
conductor. Meanwhile, the heat is removed only through 
the surface which is maintained at a fixed temperature. 
As a result, the remanent T(x,2tM) has a broad maxi- 
mum in the center, and the trapped B{x, 2tM) acquire a 
specific "bell" shape. This is illustrated in Fig. ^ for the 
case Bm — 6Bp. 

It is also interesting to analyze the evolution of the 
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FIG. 3: Trapped remanent flux after a magnetization pro- 
cess for different maximum applied fields Bm- (a) Curves for 
a — 0.48 and different r, i. e., for different rates, (b) Curves 
for r = 0.001 and different a, i. e., for different heat capac- 
ities. The solid curve is the same for both plots and shows 
the Bean-model result without any heating effects. Symbols 
show results of numerical calculations which reproduce well 
our analytical curves except for very large Bm- 

temperature and flux density during the whole magneti- 
zation process. A time sequence of curves shovifing this 
is plotted in Fig. ^, and corresponds to the intermedi- 
ate stages leading to the remanent state seen in Fig. |l| 
for T — 0.001. One can see that the temperature growth 
starts slowly and then accelerates. When the applied field 
reaches maximum (at t = Im), the T profile has already 
acquired its characteristic shape and changes little dur- 
ing the subsequent field decrease. The evolution of the 
B profile looks similar to the standard Bean-model be- 
havior. The full penetration is reached approximately at 
Ba = Bp (at t = Im/'^), and then the B{x,t) shifts up- 
ward almost uniformly until the field starts to decrease 
and fiux exits from the surface. The fiux front position 
(seen as a ridge in the 3D plot) is shifting towards the 
center, and in the final remanent state the B profile ac- 
quires its triangular shape. 

Our approach, used so far for a constant ramp rate of 
the field, is easily generalized to having different sweep 
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FIG. 4: Temperature and flux density profiles in the rema- 
nent state after application of the maximum field Bm = Bp 
(A), Bm = 1.3Bp (B), and Bm = 1.6Bp (C). The profiles 
correspond to the marked points (A,B,C) on the lowest curve 
of Fig. ^(b). The profile (B) corresponds to the maximum 
trapped flux. 



rates on the ascending and descending field branch. It 
turns out that for a given magnetization time it is some- 
what beneficial to have a faster field increase followed by 
a slower decrease to zero. In particular, for the extreme 
case of an instantaneous field increase and a descent last- 
ing 8 ms we find that one traps 8% more flux than for the 
symmetric 4-1-4 ms field pulse (for the parameters used 
in the paper). The physical reason for this is the follow- 
ing. When a larger amount of heat is dissipated in the 
beginning of the PFM, the heat has more time to flow 
out of the sample. As a result, the temperature in the 
remanent state becomes slightly lower, which results in 
a larger trapped flux. We omit to show detailed profiles, 
since these graphs are very similar to the ones already 
presented in this paper. 

The critical state model used in our analysis neglects 
the viscous force acting on fiux lines. This force can be 
an important ingredient that detemiines the B{x,t) dis- 
tributions during very fast PFM.I3 We stress, however, 
that the heating takes place independent of whether the 
flux motion is viscous or not. Our results clearly demon- 
strate that the heating produced within the critical-state 
approach is sufficient to account for the suppression of 
the trapped flux during the PFM process. 
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FIG. 5: Temperature and flux density profiles in the remanent 
state after applying magnetic field Bm ~ 6Bp. The curves 
are obtained for the field-pulse duration: 23.2 s (solid), 6.3 s 
(dashed), and 2.3 s (dotted), which correspond to r=3.2, 0.9 
and 0.32, respectively. 

IV. CONCLUSION 

The temperature and field distributions in a bulk su- 
perconductor during a PFM process are calculated ana- 
lytically within the critical state model and with account 
of heating due to flux motion. The remanent trapped flux 
<I> is smaller for a larger PFM rate, and for smaller heat 
capacity. For a given duration of the activation pulse the 
$ reaches maximum for some optimal maximum field, 
which is always smaller than twice the penetration field. 
Suprizingly, the remanent flux distribution for optimal 
field is not monotonous, but the overall peak has a small 
dip in the center. The strongest temperature rise is usu- 
ally found close to the surface. The trapped flux can be 
enhanced without changing the total PFM time if the 
field ascent is made faster than the descent to zero. 
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FIG. 6: Time evolution of temperature and flux density distr- 
butions during the full magnetization process for the maximal 
fleld Bm = 2Bp, and duration 7.3 ms (r=0.001). 



mal diffusion time), b = 7r/(2aBp), and tM is the time 
when the applied field reaches maximum, tM = Bm/Ra- 
The temperature distribution during the magnetization 
process is given by 



APPENDIX A 

Here we list analytical expressions for the temperature 
distribution for different stages of the magnetization pro- 
cess. They are obtained by substituting the heat release 
rate given by Eqs.(p|)-(pl]) and the Green function (^) 
into Eq. (^ . The expressions are given for a general case 
of having different field ramp rates Ra and Rd on the 
ascending and descending branch, respectively. The fol- 
lowing notations are used: a — tt^k/ACw'^ (inverse ther- 



T[x, t) = To H 7 > ^ Sn(t) sm — cos — — 

n=l 

(Al) 



The dimensionless coefficient S'„ are given by one of the 
expressions below, and Tab. 1 explains which expression 
should be used for every particular case. It is convenient 
to define the following functions corresponding to incom- 
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Bm < Bp 


t < tM 
tM < t < tM + Bm/Rd 


(A4) 
(A5) 


Bp < Bm <2Bp 


t < Bp/Ra 

Bp/Ra < t <tM 

tM < t < tM + Bm/Rd 


(A4) 
(A6) 
(A7) 


Bm > 2Bp 


t < Bp/Ra 
Bp/Ra < t < tM 

tM < t < tM + 'iBp/Ro 

tM + 2Bp/Rd < t <tM+ Bm/Rd 


(A4) 
(A6) 
(A7) 
(A8) 



plete and complete penetration processes: 



{anH + e""' * - l)hR 



cos{abRnt) — e 




bR/n 


— sm{abRnt) 




n2 + {hRf 





(A2) 



TABLE I: The last column shows which formula for coeffi- 



cients Sn in Eq. ( |Af| ) should be used for given Bm and t. 



Then, 



n 3 10 2 
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an^{t~Bp/R) 



S^^lt, Ra) 



Sn = Sf°™{t,RA) 

e 



e-"''(*-*^-')5r(tM, + ^ S'^'^t - tM,RD/2) 



Ka 

Sn = e-'^"'(*-*-)5™'"(tM, + ^ S';r{t - tM^Ro/^) 

Ha 



(A3) 

(A4) 
(A5) 
(A6) 
(A7) 

(A8) 
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